KIAS-P01051 
hep-ph/0111254 

Model Independent Information On Solar Neutrino Oscillations 

M.V. GarzelliQ 

Dipartimento di Fisica Teorica, Universita di Torino, 
and INFN, Sezione di Torino, 
Via P. Giuria 1, 1-10125 Torino, Italy 

C. GiuntiQ 

INFN, Sezione di Torino, 
and Dipartimento di Fisica Teorica, Universita di Torino, 
Via P. Giuria 1, 1-10125 Torino, Italy 
and 

School of Physics, Korea Institute for Advanced Study, 
Seoul 130-012, Korea 
(Dated: 13 December 2001) 

We present the results of a Bayesian analysis of solar neutrino data in terms of Ve Vti.r oscil- 
lations, independent from the Standard Solar Model predictions for the solar neutrino fluxes. We 
show that such a model independent analysis allows to constraint the values of the neutrino mixing 
parameters in limited regions around the usual SMA, LMA, LOW and VO regions. Furthermore, 
there is a strong indication in favor of large neutrino mixing and large values of Am'^ (LMA region). 
We calculate also the allowed ranges of the neutrino fluxes and we show that they are in good 
agreement with the Standard Solar Model prediction. In particular, the ratio of the *B flux with its 
Standard Solar Model prediction is constrained in the interval [0.45, 1.42] with 99.73% probability. 
Finally, we show that the hypothesis of no neutrino oscillations is strongly disfavored in a model 
independent way with respect to the hypothesis of neutrino oscillations. 
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I. INTRODUCTION 



An important model independent evidence in favor of transitions of solar electron neutrinos in other active states 
has been recently obtained from the comparison of the results of the Super-Kamiokande and SNO Q solar neutrino 
experiments 0, pi ^, ^ . This evidence is supported ||^ also by the results of the Homestake , GALLEX Q , SAGE 
1^ and GNO px|j experiments, that measured different suppressions of the solar electron neutrino flux with respect 
to the Standard Solar Model (SSM) prediction. 

The simplest and most natural mechanism that explains all solar neutrino data is v,, — > u^^t oscillations, due to 
neutrino masses and mixing (see [|l2[ [l5|). Several groups have recently analyzed the data of solar neutrino 

experiments in terms of neutrino oscillations |g ^ |2^, ||, |2|, |2|, ||, ||, |2^, H, |9). The authors 

of Refs. I, m, ^ 1^, H, H, !§ found the allowed regions in the space of the neutrino mixing parameters 
assuming that the neutrino fluxe s^ p roduced in the core of the sun are given by the BP2000 Standard Solar Model . 
The authors of Refs. jl^, HJ, H^, ^ performed analyses with the neutrino flux considered as a free parameter 
to be determined from the data (in Ref. |^ also the ''Be neutrino flux has been considered as a free parameter, 
whereas the authors of Ref. |Q considered also the hep neutrino flux as a free parameter). 

In this paper we calculate model independent allowed regions in the space of the neutrino mixing parameters. In 
the following, by "model independent" we mean independent from the solar model. In other words, we do not assume 
that the neutrino fluxes produced in the core of the sun are given by a Standard Solar Model. 
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We consider the simplest case of mixing of two neutrinos j/g and l'a^ where Va is a hnear combination of and v-r'. 

Ve = cos "i? i^i + sin -d V2 , 

Va = — sin dvi^ cos "dv^- (1) 

Here, "i? is the mixing angle and vi, vi are massive neutrinos with masses mi, mi, respectively. Neutrino oscillations 
depend on the mixing angle ■& and the squared- mass difference Am^ = — (see [ p^ |l3| , |l5|). 

Let us notice, that the treatment of solar Vf, — > Va transitions in the framework of the simple two-neutrino mixing 
scheme in Eq. (|l|) is valid with good approximation also in the case of mixing of three neutrinos (I'e, Vy. v^) with small 
C/e3 as indicated by the results of the CHOOZ long-baseline v,, disappearance experiment |Q. Furthermore, 

the treatment is also valid in 3-1-1 four- neutrino mixing schemes in which the sterile neutrino is practically decoupled 
from the active ones (1 — |C/s4p ^ 1) and J7e3 is small, as indicated |Q by CHOOZ data. 

Our statistical analysis of the data is Baycsian. We have already explained in Rcf. |^ the advantages of using 
Bayesian Probability Theory in the analysis of solar neutrino data^ (for a general introduction to Baycsian Probability 
Theory see |3 7 [ |3^ , pO| |4l[|). Another model independent analysis of solar neutrino data was performed several 
years ago [[42[ using a different statistical approach. Apart from the differences due to the statistical methods, the 
main difference with the analysis performed in Ref. is the impressive increase of the quality and quantity of data, 
that, as we will show, allow now to obtain much more stringent model independent information. 

As discussed in Ref. , the Bayesian analysis of solar neutrino data assuming the neutrino fluxes predicted by the 
Standard Solar Model takes into account the uncertainty of the fluxes in the covariance matrix of the least-squares 
function of the data, assuming a multi-normal distribution for the errors of the fluxes. This approach cannot be 
extended to the model independent analysis of solar neutrino data, in which the fluxes must be considered as free 
parameters that may be determined from the data in a similar way as the neutrino mixing parameters. Therefore, we 
have developed a new approach in which the neutrino fluxes are considered as unknown quantities with a flat prior 
probability distribution. The model independent information on neutrino mixing is calculated through marginalization 
of the posterior probability distribution with respect to the neutrino fluxes. We also derive information on each 
neutrino flux through marginalization of the posterior probability distribution with respect to the other neutrino 
fluxes and the oscillation parameters. 

In the framework under consideration, the theoretical calculation of the rates measured in solar neutrino experiments 
depends on the neutrino mixing parameters tan^ d and Am^ , on the neutrino fluxes produced in the core of the sun by 
the thermonuclear reaction of the pp and CNO cycles, and on the energy-averaged cross sections of neutrino interaction 
in the detectors. In order to have a uniform treatment of these uncertain theoretical quantities, we consider all of 
them as parameters that have some uncertainty, which is taken into account in the prior probability distribution. This 
approach is different from the standard one (see ) , in which the uncertainties of the fluxes and cross sections are 
taken into account in the covariance matrix that determines the likelihood function (sampling distribution). However, 
we think that the new approach is equally, if not more, reasonable than the standard one, because the values of the 
fluxes and cross sections are part of our prior knowledge. In any case, the two approaches should give approximately 
the same result in the model dependent analysis in which the neutrino fluxes are assumed to be those predicted by 



the Standard Solar Model. We check this consistency in Section IV 



In the following it is convenient to group the solar neutrino fluxes in the array 

$ = ($1,... (2) 

The neutrino flux <I>i is produced in the i**^ thermonuclear reaction, with z = 1, . . . , 8 for the pp, pep, /lep, ^Be, ^B, 
i3]\j^ 15q^ i7p thermonuclear reactions, respectively. 

In both the model dependent and independent analyses we perform a Rates Analysis that takes into account the 
total rates measured in the Homestake, GALLEX-I-GNO-I-SAGE, SNO, Super-Kamiokande experiments, and a Global 
Analysis in which the Super-Kamiokande total rate is replaced by the 38 bins of the Super-Kamiokande day and night 
electron energy spectra. Hence, calling A^s the number of solar data points, we have A'^s — 4 in the Rates Analysis 
and A'g = 41 in the Global Analysis. We call the energy-averaged cross section of detection of the i*'' neutrino 
flux for the j"^ solar data point, with j = 1, . . . , A^s- We group these cross sections in the 8x As matrix 

(3) 




Bayesian Probability Theory has been appUed to the analysis of solar neutrino data also in Refs. |l9| p^ ]. 
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Bayesian Probability Theory allows to calculate the posterior probability distribution p(tan^??, Arn^, <J>, C|D, I) for 
the neutrino mixing parameters, the neutrino fluxes and the detection cross sections through Bayes Theorem: 

,(tan^., A^^<i>,C|D,I) ^ MD| tan^^, C, I) ,(tan^^, A^-|I) ,(cl>|I),(C|I) 

p(D|I) 

where D represents the data and I represents all the background information and assumptions on solar physics, 
neutrino physics, etc. The sampling distribution p(D| tan^i?, Am^, $, C, I) is also known as likelihood function. The 
prior distributions p{tan^'&, Am^ |I), J3(<1>|I) andp(C|I) quantify our prior knowledge, or lack of knowledge, on the values 
of the mixing parameters, neutrino fluxes and detection cross sections, respectively. The probability distribution p(D| I) 
is known as global likelihood and acts as a normalization constant. 

The prior probability distribution of mixing parameters, p(tan^?9, Am^|I), must be chosen in order to quantify 
appropriately the prior lack of knowledge on the values of the parameters tan^i? and Am^. Assuming ly^ —>■ v^^t- 
oscillations (included in I), we know that solar neutrino data are sensitive to several different order of magnitude of 
tan^t? and Am^ , through vacuum oscillations or resonant MSW transitions . Therefore, the most reasonable 

non-informative prior probability distribution function, that we will use in the following, is a flat distribution in the 
log(tan^t9)-log(ATO^) plane. Indeed, the authors of Ref. have chosen, independently, the same prior, which 
is consistent with the calculation of the credible regions for the neutrino oscillation parameters in the log(tan^i?)- 
log(ATO^) plane. 

Using a flat prior probability distribution in the log(tan^??)-log(ATO^) plane, Eq. becomes 

rt. 2, . 2 ^ PIT. .^ p(D|tan2^,Am^<l>,C,I)p(<i>|I)p(C|I) 

p^ran ism , i) j^i^^^^^^2^^ dlog(Am2) p(D| tan^,?, Am^ $, C, I) p($|I) p{C\l) ' 

where we have expressed the global likelihood p(D|I) as the appropriate normalization constant. Integral probabilities 
must be calculated integrating p{ta,ia^-&, Am?, $, C|D, I) over dlog(tan2?9) dlog(A77i2) d<i> dC. 

The prior probability distribution of cross sections, p(C|I), represents our knowledge of the cross sections, which is 
reasonabl y ac curate. Therefore, we take this prior probability distribution as a multinormal distribution as explained 



in Section III 



As explained in Subsection IV A , in the Standard Solar Model dependent analysis the prior probability distribution 



of neutrino fluxes, p($|I), is a multinormal distribution. In the model independent analysis we consider a flat prior 



probability distribution of neutrino fluxes, as explained in Subsection VA. 

The posterior distribution for the oscillation parameters tan^ t?, Arrr \s obtained through marginalization of the 
posterior distribution (^) with respect to the neutrino fluxes and cross sections: 

p{t&n^d,Am^\D,l) ^ j j AC p{t&n^ d , Am^ , ^ , C\D ,1) . (6) 

Using this posterior distribution, we calculate Bayesian credible regions^ in the plane of the oscillation parameters 
tan^i? and Am^ 

In both the Rates Analysis and the Global Analysis we consider also the data of the CHOOZ experiment, that are 
important because they exclude large mixing for Am^ > 10^'^ eV^ p^ . 

For convenience, in the discussion we use the standard SMA, LMA, LOW and VO names for the regions in the 
tan^t^-Am^ plane. We consider these regions as enlarged with respect to the usual ones (see [^), because the model 
independent credible regions that we obtain are slightly larger than the corresponding model dependent ones: Small 
Mixing Angle (SMA) for 

tan^T? < 10"^ , 10"*^ eV^ < Am^ < 10"^ eV^ , (7) 

Large Mixing Angle (LMA) for 

10^2 < ia.T?d < 10^ 3 X 10"^ eV^ < Am^ < 10"^ eV^ , (8) 

LOW for 

10^2 < tan^^? < 10^ , 10~^ eV^ < Am^ < 3 x 10~^ eV^ , (9) 



Credible regions, also known as highest posterior density regions, contain a specified fraction of the posterior probabilit y a nd all values 
of the parameters inside the credible regions have higher probability density than those outside (see, for example, Rcf. |38l). 
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Vacuum Oscillation (VO) for 

10"^ < ta.T?d < 10^ , Am^ < 10"® eV^ . (10) 

The plan of this paper is as follows. In Section ^ we present the sampling distributions in the two types of analysis 
that we perform: the Rates Analysis of Homestake, GALLEX+GNO+SAGE, SNO and Super-Kamiokande total rates 
and the Global Analysis of Homestake, GALLEX+GNO+SAGE and SNO rates and Super-Kamiokande day and night 



spectra. In Section III we present the prior distributions for the cross sections in the Rates and Global Analyses. In 
Section [V we perform a SSM model dependent analysis in order to check the validity of our approach by comparing 
the results with those obtained with a standard Bayesian analysis in Ref . |^ . In Section |^ we present and discuss 
the results of the model independent analysis of solar neutrino data. In the concluding Section W\ we summarize our 
results. 

II. SAMPLING DISTRIBUTION 

F ollow ing the tradition, we perform two types of analysis, with different data sets. In the Rates Analysis (Subsec- 



tion II A ) we consider the total rates measured by solar neutrino exp eriments (Homestake, GALLEX+GNO+SAGE, 
SNO, Super-Kamiokande). In the Global Analysis (Subsection H B ) we consider the total rates of the Homestake, 
GALLEX+GNO+SAGE, SNO experiments and the Super-Kamiokande day and night electron energy spectra. Both 
analyses take into account also the data of the CHOOZ experiment ||33|| . 

Since the solar and CHOOZ data are independent, the sampling probability distribution can be written as 

p(D| iar?d, Am^ $, C, I) = p(Ds| tan^i?, Am^ C, I) p(Dc| tan^i?, Am^ I) , (11) 

where Ds represents the solar data and Dc represents the positron spectra measured in the CHOOZ experiment, that 
are obviously independent from the solar neutrino fluxes and detection cross sections. 

Assuming a normal distribution of experimental errors, the sampling probability distribution of CHOOZ data is 
given by 

p(Dc|tan2z?,Am2,I) = /dac ^"^^^ , (12) 

where we have marginalized over the nuisance parameter ac, that is the absolute normalization constant of CHOOZ 
positron spectra ISSj. Here Nc = 14 is the number of CHOOZ data points, is the CHOOZ least-squares function 
and Vc is the corresponding covariance matrix of uncertainties. 

The CHOOZ least-squares function is calculated as in the analysis A of Ref. |3^, with the following approx- 
imations. Since we do not know the antineutrino spectrum, the spatial distribution functions of the reactor cores 
and detector and the detector response function linking the real and visible positron energies, for each energy bin we 
calculated the oscillation probability at the average energy of the bin and at the average distance of the detector from 



each of the two reactors 1 47 . This approximation is quite good because we are interested in small values of Am^ , for 



which the energy and distance dependence of the survival probability of the z^e's in the CHOOZ experiment is very 



weak. We calculate X^ as in Eq. (13) of Ref. |33 , with the only difference that we neglect the energy-scale calibration 



factor, whose small uncertainty (1.1%) is practically negligible 



«c ) (^c ^).U2 {iC - «c ) + (^) ' (13) 



where = 2.7 x 10"^ is the uncertainty of the absolute normalization constant ac |3^- We calculate the CHOOZ 
covariance matrix Vc as described in Eq. (12) of Ref. [Q. The only missing information in Ref. |Q is the value of 
the systematic uncertainties of the positron energy bins, for which only the values for positron energy 2 and 6 MeV 
are given. For the other bins we take systematic uncertainties interpolated linearly between these two values. 

A. Rates Analysis 

Assuming a normal distribution of experimental errors, the sampling probability distribution of solar data in the 
Rates Analysis is given by 

p(Ds|tan2^,Am2,<D,C,I) = — 4-1^=, (14) 
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where Ns = 4 is the number of solar data points. X| is the solar least-squares function and Vs is the corresponding 
covariance matrix. 

The solar least squares function Xg is given by 

^l=E ^ ,2 ^ ^ (15) 

where Rj"^'^ is the event rate measured in the j"^ experiment and i?^*''^ is the corresponding theoretical event rate, that 
depends on Am^, tan^O, the neutrino fluxes and the neutrino interaction cross sections in the detectors. The index 
j = 1, . . . ,4 indicate the four types of solar neutrino experiments listed in Table | together with the corresponding 
event rates and experimental uncertainties CTj, calculated by adding in quadrature the statistical and systematic 
uncertainties for each experiment. 

Since the uncertainties of the rates of different experiments are uncorrelated, the determinant |Vs| of the solar 
covariance matrix in Eq. (|l4|) is simply given by 



Ns 

i^si=n'^i- (16) 



The theoretical event rate R^ can be written as 



^th)^^^(th)^ (17) 



where 



i?! = $z a, {Am' , tan' 9) , (18) 



is the theoretical event rate in the j*^ experiment due to the neutrino flux <i>i produced in the i^^ thermonuclear reaction 
of the pp and CNO cycles in the sun (i = 1, . . . ,8 for the pp, pep, hep, ^Be, ®B, ^^N, ^^O, ^^F thermonuclear reactions, 
respectively). is the corresponding energy-averaged cross section and (Am^, tan^ 0) is the corresponding 

averaged survival probability of solar i^e's, that depends on Ato^ and tan^ 6' (in the case of the Super-Kamiokande 
experiment, j = 4, also the averaged i^e Va transition probability must be properly taken into account). 

For the calculation of the probabilities Pij(Am^, tan^ 0) we have used the tables of neutrino spectra, solar den sity 
and radiochemical detector cross sections available in Bahcall's web pages W8 , BP2000 Standard Solar Model [hi . 
For the calculation of the theoretical rate of the SNO experiment we used the charged-current cross section given 
in Refs. |5^, 0. We did not consider the contribution of radiative corrections discussed in Refs. |5|, which 
increase the neutrino-deuteron cross section by a few percents, that is smaller than the experimental uncertainty of 
the SNO event rate ||, ||. 

The probability of neutrino oscillations is calculated with an unified approach that allows to pass continuously from 
the vacuum oscillation regime to MSW transitions Q, ^ through the quasi- vacuum regime j55[ |5^ , using the quasi- 
vacuum analytical prescription given in Ref. (see also Refs. ^), the usual prescription for the MSW survival 
probability (see j5^, Q) and the level crossing probability appropriate for an exponential density profile [6^ . We 
calculate the regeneration in the Earth using a two-step model of the Earth density profile |]63| , pi] , |65| , |66[ |6711 , that is 
known to produce results that do not differ appreciably from those obtained with a less approximate density profile. 



B. Global Analysis 

In the Global Analysis, instead of the total rate of the Super-Kamiokande experiment we consider the data on the 
Super-Kamiokande day and night electron energy spectra presented in Ref. , that contain information on the total 
rate plus the shape of the energy spectrum. Assuming a normal distribution of experimental errors, the sampling 
probability distribution of solar data in the Global Analysis is given by 

p(Ds|tan^z^,Am^<I>,C,I)= ^ ^ ""Z'^fnT-^ ^ (1^) 

(27r)^s/2y|V^ (27r)^sK/2^|ygj^| 
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where the first factor takes into account the rates measured in the Homestake, GALLEX+GNO+SAGE and SNO 
experiments. Therefore, the first factor in Eq. is calculated as in the Rates Analysis, with the only difference 
that iVs = 3. 

The second factor in Eq. ( |l9| ) takes into account the Super-Kamiokande data on the day and night electron energy 
spectra. 

The Super-Kamiokande least-squares function Xg-^ is given by 



38 

^SK = 



with the indexes 6i, 62 run from 1 to 38 for the bins of the Super-Kamiokande day and night electron energy 
spectra given in Table III of Ref. ||6^. The indexes 61,62 = 1, . . . ,19 refer to the day spectrum and the indexes 
61, 62 = 20, ... ,38 refer to the night spectrum. The covariance matrix Vsk is written as 

The statistical uncertainties ((T^''*^-')b are given in the third (for b = 1,... ,19, day spectrum) and fourth (for b — 
20, . . . ,38, night spectrum) columns of Table III in Ref. |6^. The experimental uncorrelated systematic uncertainties 
(o'u'Sc ^)(, (for 6 = 1,... , 38) are given by '^I^unc ^^^^ listed in the third (for 6 = 1, . . . , 19, day spectrum) 

and fourth (for 6 = 20, . . . , 38, night spectrum) columns of Table III in Ref. [^. The last term in Eq. ( |2l| ) is due to 
the assumption of a full correlation of the systematic uncertainties of the day and night bins with the same energy. 
The values of <5^™^^ = ^I+ig unc for 6 = 1, . . . ,19 are listed in the sixth column of Table III in Ref. (we took the 
bigger between the positive and negative values). 

The experimental correlated systematic uncertainties are given by (CTcor''')f, = Rb^^^ ^i^col- Unfortunately, the 
relative correlated systematic uncertainties presented by the Super-Kamiokande collaboration in the fifth column of 
Table III in Ref. |68| contain the contribution from the energy-dependent part of the theoretical uncertainty on the 
®B neutrino spectrum given in Ref. |^^, that is taken into account by the prior probability distribution discussed in 
Subsection [II B , We list in Table || the values of relative theoretical uncertainty i5^*cor f'^'' ^^'^ rates in the Super- 



Kamiokande spectral bins (6 = 1, . . . ,38) due to the theoretical uncertainty on the B neutrino spectrum extracted 
from Fig. 5 of Ref. [ |69) . In order to avoid double-counting, we subtracted quadratically the energy-dependent part 
of {i.e. (5^*''^^, — ^i^coiJ where 6 = 1 is the index of the first Super-Kamiokande energy bin) from the relative 

correlated systematic uncertainties in the fifth column of Table III in Ref. ||6^ (for which we took the bigger between 
the positive and negative values). The resulting values of Sf^ are listed in Table ||. 



III. CROSS SECTIONS 



In this section we discuss the prior distribution p(C|I) for the energy-averaged cross sections of neutrino interaction 
in the detectors. In the case of the Super-Kamiokande experiment there is practically no uncertainty on the value 
of the neutrino-electron elastic cross section, but there is an uncertainty of the shape of the *B neutrino spectrum 
p9| , |70|] . We take this uncertainty into account both in the Rates Analysis and Global Analysis, as described in the 
following subsections. 



A. Rates Analysis 

In the Rates Analysis the number of experimental data points is Ns — 4 (see Table |), leading to a 8 x 4 matrix C 
of energy averaged cross sections. 

Assuming a complete correlation of the errors of the energy-averaged cross sections of different neutrino fluxes 
in each experiment ]7l[ |, the cross section of neutrino detection in the j**^ experiment averaged over the energy 
spectrum of the i^^ neutrino flux can be written as 

C., =(l + e,Alna,)CO , (22) 

where C,° is the standard averaged cross section, AlnC^ is the relative uncertainty of C^j, and is an unknown 
parameter that quantifies the deviation of the averaged cross section from its standard value. Assuming a normal 
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distribution of errors, for each we adopt the prior distribution 

*'I" = W- '''' 

The prior distributions for different 's are independent because the theoretical errors of the energy-averaged cross 
sections of neutrino interaction in different experiments are uncorrelated. The prior distribution p(C|l) of energy 
averaged cross section is obtained from Eqs. ( p^ ) and (|2^). In practice, since all quantities that we calculate are 
integrated over dC using a Monte Carlo, the prior distribution p(C|l) is generated by the Monte Carlo from the prior 
distribution (|2|) of , using Eq. (|2|). 

This procedure may seem arbitrarily complicated, but one must notice that it is not possible to write a multi-normal 
distribution for fully correlated quantities as Ctj for i = 1, . . . ,8, because the covariance matrix would be singular. 
Indeed, complete correlation means that the errors are not independent: the deviation of one Cij with respect to C'j'j 
determines the deviations of all the others Ci'j from C^,j, with i' ^ i, as quantified in Eq. (^). 

The values of C^j that we use for the Chlorine (j = 1) and Gallium {j = 2) experiments are given in Table V of 
Ref. ||3^. For the SNO experiment (j = 3) with kinetic energy threshold T^^° = 6.75 MeV for the detected electrons, 
we calculated 

Cfg = 1.82 X lO^^cm^s, C^g = 1.98 x 10"^ cm^ s , (24) 

and all the other Cf^ are zero. For the Super-Kamiokande experiment {j — 4) with kinetic energy threshold — 
4.50 MeV for the detected electrons, we calculated 

C^4 = 1.04 X lO^^cm^s, = 1.98 X IQ-^cm^s, (25) 

and all the other Cf^ are zero. 

The values of the relative uncertainties AlnC^ = ACij/Cij for j = 1,2 ('^^Cl, ^^Ga experiments) are given in 
Table VI of Ref. ||. For the SNO experiment {j = 3), AlnCgg = AlnCga = 3.0 x lO^^ |, |9l ^ and all the other 
AlnCis are zero. For the Super-Kamiokande experiment {j = 4) we take AlnC54 = 2.0 x 10^^^ in order to take into 
account the uncertainty of the shape of the ®B neutrino spectrum calculated in Ref. All the other AlnCi4 are 
zero. 



B. Global Analysis 

In the Global Analysis the prior distribution of the averaged cross sections in the Chlorine (j = 1), Gallium (j = 2) 
and SNO {j = 3) experiments is the same as in the Rates Analysis. For the Super-Kamiokande experiment instead 
of the total rate, we consider the data on the 38 bins of the Super-Kamiokande day and night electron energy spectra 
(see Table Hence, the number of experimental data points in the Global Analysis is Ng = 41, leading to a 8 x 41 
matrix C of energy averaged cross sections. 

The energy-averaged cross sections Cij for j < 3 are given by Eq. (|2^), with the distribution of the normalized 
deviations £,j given in Eq. (p3|). 

The energy-averaged cross sections C^j, with j = 4, . . . , 41, of neutrino detection for the Super-Kamiokande day 
and night spectral bins that take into account the fully correlated uncertainty due to the uncertainty of the shape of 
the neutrino spectrum |^, |7^ are given by 

^5(6+3) = {l+U C) QV+3) (6 = 1, . . . , 38) , (26) 

with the normahzed deviation ^4 distributed as p(C4|I) given in Eq. ( p3| ) with j = 4. In Eq. (|2^), C^i^f^^^^^ is the 

standard energy-averaged cross section of ®B neutrino detection in the 5'^ Super-Kamiokande bin given in Table 0, 
together with C'^j-^^g-j , relative to the hep neutrino flux (we take 6*3(^+3) = Cg^^^g-) and all the other 6*^(^+3) = 0). 

The values of relative theoretical uncertainty S^*"^^^ given in Fig. 5 of Ref. j6^ for the rates in the Super-Kamiokande 
spectral bins (6=1,... , 38) due to the theoretical uncertainty on the ^B neutrino spectrum are listed in Table ||. 



IV. MODEL DEPENDENT ANALYSIS 



In this Section we perform a Bayesian model dependent analysis of solar neutrino data assuming the solar neutrino 
fluxes predicted by the BP2000 Standard Solar Model [jll]. In this analysis we take into account of all the theoretical 
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uncertainties through the prior distributions. We check the vahdity of this approach comparing the results with those 
obtained with a standard Bayesian analysis in Ref . |^ . 

The sampling distributions and the prior dist ributi ons of detection cross sections are described, respectively, in 
Sections || and II] . In the following Subs ection IV A we describe the SSM prior distribution of fluxes used in the 
model dependent analysis. In Section tVB we discuss our results for the calculation of the Bayesian credible regions 
in the plane of the oscillation parameters tan^i^ and Am^, to be compared with the corresponding ones presented in 
Ref. [p5l. 



A. SSM Prior Distribution of Fluxes 



In the model dependent analysis we assume a multi-normal prior distribution for the neutrino fluxes centered on 
the BP2000 SSM values of the neutrino fluxes, ^p'^ with the standard correlated uncertainties (see ^): 



p($|I) = N cxp 



-I E i'^n - '^^n (v^'')n 



■12 (*»2 



SSM\ 



ii ,12 = 1 



(27) 



where TV is a normalization constant and 6'($i) — 1 for $i > and 9{^i) — for $i < 0. V$ is the covariance matrix 
of the fluxes, given by U% 



12 



(28) 



k=l 



Here X}~, with /c = 1, . . . , 12, are the twelve input astrophysical parameters in the SSM (^i^i, S's^s, >S'3_4, <S'i,i4, S'l,?, 
Luminosity, Z/ X , Age, Opacity, Diffusion, Ctbc, •S'l.ie); whose relative uncertainties AlnX^, determine the correlated 
uncertainties of the neutrino fluxes (I)?^'^ through the logarithmic derivatives 



Oiik 



ain$fM 
ainXfe 



(29) 



We adopt the values of a^jt and AlnXj, given in Ref. [Q, except for AlnX7 (relative to Z/X), whose value has been 
updated in the BP2000 SSM Q from 0.033 to 0.061, and ai,i2 = i5i8, AlnXia = 0.181, that have been introduced 
for the first time in the BP2000 SSM X12 is the astrophysical factor S'i,i6 for the reaction i^0(p,7)"F that 

determines the small ^^F neutrino flux (i = 8). 



B. Credible Regions 

The credible regions with 90%, 95% and 99% posterior probability calculated through Eq. (^ are shown in Fig. |l| 
for the Rates Analysis and in Fig. || for the Global Analysis. One can see that they are almost identical to the 
corresponding regions in Figs. 3 and 4 of Ref. ||2^, obtained with the standard method, in which the uncertainties of 
the fluxes and cross sections are taken into account in the covariance matrix that determines the likelihood function. 



The integral probabilities of the LMA, VO, SMA and LOW regions presented in Table |I| almost coincide with 
those obtained in Ref. As noted in Ref. ||2^, these probabilities favor the LMA solution of the solar neutrino 

problem. 



Table Inil shows also the values of the integral probabilities of large and small values of tan^t? and Am^, that favor 
large mixing and large values of Am^. 

In conclusion of this section, we have shown that in the model dependent analysis of solar neutrino data the new 
approach that we follow here, in which all the theoretical uncertainties are taken into accounts in the priors, gives 



approximately the same results as the standard one followed in Ref. |25 , in which the theoretical uncertainties are 
taken into account in the covariance matrix that determines the likelihood function. Therefore, we conclude that the 
new approach is compatible with the standard one and we can apply it to the model independent analysis of solar 
neutrino data. 

Let us notice, however, that the new approach is not convenient for the model dependent analysis of solar neutrino 
data, because it needs the integrations over the neutrino fluxes and cross sections in Eq. (||). 
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V. MODEL INDEPENDENT ANALYSIS 



This is the main Section of the paper, in which we present the resuhs of the model independent analysis of solar 
neutrino data assuming a flat prior probability distribution of neutrino fluxes, that is explained in the following 
Subsection VA. 



Am . In Subsection 
posterior distribution 



In Subsection VB we present the model independent allowed ranges for the neutrino mixing parameters tan i?, 

2 



Finally 

oscillations and neutrino oscillations 



V C we discuss the information on the solar neutrino fluxes obtained from the model independent 
in Subsection VD we compare the probabilities of the hypotheses of no neutrino 



A. Model Independent Prior Distribution of Fluxes 

The Standard Solar Model has reached a reasonable degree of credibility in recent years as a consequence of its 
very good agreement with helioseismological measurements Therefore, we think that the true values of the 

neutrino fluxes cannot be too different from those of the SSM. 

In our model independent analysis we assume a flat prior distribution from zero to twice the BP2000 SSM value 
im for all the neutrino fluxes except hep and pp. 



7j), the value of the 
le Super-Kamiokande 



Since the astrophysical ^-factor of the hep reaction is very difficult to calculate (see Ref. [ 
hep flux is the most uncertain one and no uncertainty is given in the BP2000 SSM. However, t^ 

Collaboration measured Q a 90% CL upper limit for the hep neutrino flux of 4.3 times the BP2000 SSM prediction 
Hence, in the model independent analysis we take for the hep flux a flat prior distribution from zero to ten times the 
BP2000 SSM value 

The prior distribution of the pp flux is obtained from the prior distributions of the other fluxes and the luminosity 



constraint (see Ref. [75 



^i^Q, (30) 

i 

where a; is the average thermal energy released together with a neutrino from the source i given in Ref. ||75|| and 
Kq = Cq/AttR^ = 8.527 x 10" MeVcm-^ s~i is the solar constant (£© = 2.398 x 10^^ MeVs~^ is the luminosity of 
the sun and R = 1.496 x 10^^ cm is the sun-earth distance). The choice of the pp flux as the dependent flux to be 
determined by the luminosity constraint (|3^) is mandatory, because the pp flux, being the largest one in the SSM, is 
the only flux that can exhaust by itself the luminosity constraint without assuming an unreasonably large value. 

The resulting prior distribution of the ratio $pp/<I>pp^ is almost flat in the range [0.904, 1.094], as can be seen from 
Figs. |l5| and |l^ (thick dashed curve). The upper value of this range corresponds to the maximum pp flux allowed 



by the luminosity constraint (30). The lower value depends on the choice of the range of the prior distributions of 
the other fluxes. As shown in Figs. ^ and |l^, the prior distribution of the pp flux is almost exactly symmetric with 
respect to $pp/$|p^ = 1 and is exactly flat in the range [0.94, 1.06] (the small fluctuations of the prior distribution 
in Figs. |l^ and |l6| is due to the fact that the prior distributions of the neutrino fluxes are generated through a Monte 
Carlo, in order to perform the integrals over the neutrino fluxes in Eq. (H)). 

In Table we list the relative prior uncertainties A$i/$p^ of the neutrino fluxes in the model independent 
analysis and we confront them with the relative uncer tainties {M^^^ / ^^^^) in the BP2000 Standard Solar 
Model and with the relative differences between the ffuxes predicted by several old and new models and the BP2000 
SSM. One can see that the relative prior uncertainties in the model independent analysis are much larger than the 
relative uncertainties in the BP2000 Standard Solar Model and they are also much larger than the relative differences 
between the fluxes predicted by other models and the BP2000 SSM. Therefore, our choice of prior distributions for 
the neutrino fluxes is indeed appropriate for a model independent analysis of solar neutrino data. 



B. Oscillation Parameters 



In this Subsection we present our results for the model independent credible regions in the tan^i9-Am^ plane 
obtained with the posterior distribution (^ marginalized with respect to the neutrino fluxes and cross sections. 

The credible regions with 90%, 95% and 99% posterior probability obtained with the Rates Analysis are shown 
in Fig. p| One can see that they are significantly larger than the model dependent credible regions in Fig. n^, but 
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still restricted in the parameter space: a large area of the parameter space shown in Fig. S is excluded^ in a model 
independent way. In particular, the absence of neutrino oscillations, corresponding to (Att? tan^ t?) < 10~-'^^eV^, is 
excluded. 

In Fig. Ijone can distinguish clearly the 90% and 95% allowed LMA, SMA, LOW and VO regions that are obviously 
larger than the corresponding regions obtained in the model dependent analysis (Fig. |^). The 99% allowed regions 
are less defined, but large values of Ato^, corresponding to the LMA and SMA regions with Am^ ^ 3 x lO^^eV^, are 
separated from small values of Am^ , corresponding to the LOW and VO regions with Am2 < 3 X lO-^eV^ The LMA 
and SMA 99% allowed regions are connected through a channel at Am^ ~ 1.3 x lO^^'eV^. The LOW 99% region is 
connected with the upper VO 99% region, but there are still two separated VO 99% regions at Am^ ~ 3 x 10^^°eV^ 
and Am2 ~ 1 x 10-i°eV^ 

All experiments are important in order to restrict the model independent allowed regions, as one can see confronting 
Fig. ^ with the four figures ^-Q, each one obtained neglecting the data of one of the four experiments in Table |. 
In particular, one can see that the Gallium data of GALLEX+GNO+SAGE and the rate measured in the Super- 
Kamiokande experiment are crucial in order to separate the 90% LMA and LOW regions, and all experiment are 
important for the separation of the 95% and 99% LMA and LOW regions. The SNO and Super-Kamiokande rates 
are important for the exclusion of a large part of the "dark-side" region (tan^ ■&>!). 

Figure ^ illustrates the quality of the fit of solar neutrino data in the Rates Analysis for some selected values of 
tan^i? and Am^. The corresponding best-fit values of the neutrino fluxes, normalized to the BP2000 SSM prediction, 
are given m Table 0. 

The solid, long-dashed, short-dashed and dotted fines in Fig. ^ correspond, respectively, to the best-fit points in 
the SMA, LMA, LOW and VO regions (we call "best-fit point" the point with highest posterior density of probability) . 
From Fig. ^A one can see that the best-fit LMA solution fits very well the data, the best-fit SMA solution fits well 
the data, but it is slightly off the SNO rate, the best-fit VO solution do not fit well the Gallium rate, and the best-fit 
LOW solution fit of the Gallium and Chlorine data is rather bad. 

The solid line in Fig. ||B corresponds to tan^i? = 1.0 x 10~^ and Am? = 3.0 x 10^^ eV^, a point near the center of 
a quasi-triangular excluded region in Figs. From Fig. one can see that in this case the energy dependence of 
the oscillation probability is very strong, typical of small mixing angle MSW transitions (see ||l^, |l^, |l^, 45 1), with 
small suppression of low-energy pp neutrinos that do not cross the resonance and a large suppression of intermediate 
and high energy neutrinos. The fit of the Super-Kamiokande rate is less disastrous than the fit of the SNO rate 
because of the high flux of f^_r due to a large *B flux at the limit of the maximum value <I>SB/$f™ = 2, as shown in 
Table 0. From Fig. one can infer that the quasi-triangular excluded region centered at tan^?9 ~ 1.0 x 10 ^ and 
Am^ ~ 3.0 X 10~^ eV is still excluded eliminating one data point, in agreement with Figs. ^-0. 

The long-dashed line in Fig. ^ corresponds to tan^i^ = 5.0 x 10~^ and Am^ = 1.0 x 10^^ eV^, a point in a 
quasi-triangular excluded region at the center of Figs. ||-^ One can see that in this region it is not possible to fit 
the Gallium, Chlorine and SNO data. The reason is a strong suppression of all energy neutrinos due to the crossing 
of the MSW resonance in the sun. Only the Super-Kamiokande rate can be fitted well taking a large ^B flux at the 
limit of the maximum value of the prior distribution, as shown in Table M. 

The short-dashed line in Fig. ||B corresponds to tan^-d = 1.00 and Am — 9.0 x 10^^^ eV^, a point near the center 
of an excluded area in the VO region in Figs. ||-^ and 0. One can see that this region is excluded by the SNO rate, 
in agreement with the fact that it is not excluded in Fig. |^, obtained with the model independent Rates Analysis of 
Gallium, Chlorine and Super-Kamiokande data only. 

The dotted line in Fig. |B corresponds to tan^i? = 1.00 and Am^ = 5.5 x lO^^^eV^ a point in the "Just So^" 
region found in Ref. in the SSM model dependent analysis of pre-SNO solar neutrino data. One can see that the 
SNO rate disfavors this region, that is indeed allowed in Fig. ||, obtained with the model independent Rates Analysis 
of Gallium, Chlorine and Super-Kamiokande data only. 

Finally, the dash-dotted line in Fig. ^ corresponds to the case of no oscillations. Obviously, as shown in Table 0, 
in this case the flux of intermediate energy neutrinos must be suppressed and the *B neutrino flux must have an 
intermediate value between those measured in the Super-Kamiokande ^ and SNO |^ experiments assuming no 
neutrino oscillations. The smaller experimental uncertainty of the Super-Kamiokande rate forces the ^B neutrino flux 
to be close to the one measured in the Super-Kamiokande experiment in the case of no-oscillations, leading to a bad 
fit of the SNO rate. This ^B neutrino flux adds to thepp flux to give a rate in Gallium experiments that is too high, 
leading to a bad fit of Gallium data, as shown in Fig. [qB. 

Let us consider now the Global Analysis. Figure pTshows the credible regions with 90%, 95% and 99% posterior 



For convenience, we call "excluded" the areas outside of the credible regions, although they are not strictly speaking excluded, since 
they have a non-zero, albeit small, posterior probability. Similarly, we sometimes call "allowed" the areas within the credible regions. 
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probability obtained with the Global Analysis. Again, the allowed regions are larger than those obtained with the 
model dependent analysis (see Fig. ||) , but they are still clearly distinguishable and a large region of the parameter 
space is excluded. One can notice that in the model independent Global Analysis there is a very small SMA credible 
region at the 99% probability level, that is not present in the model dependent analysis. 

In Fig. I the LMA regions with 90%, 95% and 99% probability are clearly separated from the other regions, whereas 
the LOW and VO credible regions are connected through the quasi- vacuum region at Am2 - 6 X 10"^ eV^ Several 
disconnected VO credible regions are present for 8 x 10~^^ eV^ < /S.m? ^ 3 x 10"^ eV^. 

Figure ^ illustrates the quality of the fit of solar neutrino data in the Global Analysis for some selected values of 
tan^-(? and Am^, with the corresponding best values of the neutrino fluxes normalized to the BP2000 SSM prediction 



given in Table VI. 

The lines in Fig. p^|A correspond to the same values of tan^i? and Am^ as in Fig. ||A, i.e. to the best- fit point in the 
SMA, LMA, LOW and VO regions obtained in the Rates Analysis. The solid, long-dashed, short-dashed and dotted 
lines in Fig. |10|B correspond, respectively, to the best-fit point in the SMA, LMA, LOW and VO regions obtained in 
the Global Analysis. The dash-dotted line in Fig. [io|B corresponds to the case of no oscillations. 

The best-fit SMA solution in the Rates Analysis (solid line in Fig. ^^|A) is disfavored by the Super-Kamiokande 
energy spectrum, because the corresponding electron- neutrino survival probability increases with energy too steeply. 
The best-fit SMA solution in the Global Analysis (solid line in Fig. [1C|B) has a smaller mixing angle, leading to 
a flatter energy dependence of the electron-neutrino survival probability, in agreement with the Super-Kamiokande 
spectrum. The survival probability of high-energy electron-neutrinos is close to one and the Super-Kamiokande data 
are fitted with a small *B neutrino flux, as shown in Table [V^ ( <I>8b/'I'sp'^ — 0.53). In this case the theoretical rate 
for the SNO experiment (j = 3) is too high, as shown in Fig. |10p3, disfavoring the SMA solution. 

Figure |l^ shows that the data are well fitted in the LMA and LOW regions, whereas the VO solution is slightly 
disfavored by the highest-energy bins of the Super-Kamiokande spectrum. It is interesting to notice the difference of 



the LMA prediction for the day and night Super-Kamiokande energy bins, that is especially evident in Fig. 10 A. This 
difference is due to the regeneration of electron neutrinos in their passage through the Earth. 

The case of no oscillations (dash-dotted line in Fig. 10|B) is compatible with the Super-Kamiokande energy spectrum 
through a small *B neutrino fiux, as shown in Table V| (<I>8b/$8™ — 0.44). However, this fiux is too high to fit the 
Gallium {j = 1), Chlorine (j = 2) and SNO [j = 3) rates. 

Let us now determine the probabilities of the SMA, LMA, LOW and VO regions, that are given by the integrals of 
the posterior density over the appropriate ranges of the parameters given in Eqs. (|7|)-([lO|): 

p(R|D,I)= / dlog(tan2i5)dlog(Am2)p(tan2,9, Am2|D,I) , (31) 

with R = SMA, LMA, LOW, VO. The values of these probabilities in both the Rates Analysis and Global Analysis 
are listed in Table III, together with those obtained in the model dependent analysis presented in Section IV (see 
also ||2^). One can see that the model independent and model dependent integral probabilities of the four regions 
are similar and the LMA region is favored in both cases, especially with the Global Analysis that, on the other hand, 
strongly disfavors the SMA region. 

Bayesian Probability Theory allows also to calculate the separate posterior probability distributions for ia.T?d and 
Am^ through the marginalizations 

p(tan2i?|D,I) = y"dlog(Am2)p(tan27?,ATO2|D,I), (32) 

p(Ato2 |D, I) = y dlog(tan2i9) p{i'AT?d, Ato^ |d, I) . (33) 

Figures |ll| and 12 show the model independent posterior distribution of tan^i? in the Rates Analysis and Global 
Analysis, respectively. One can see that in both cases there is a wide and high peak at large mixing angles (tan^?9 ^ 
0.3). In the Rates Analysis there is also a small peak at small mixing angles (tan^i? ^ 10~^). The integral probabilities 



of large and small mixing are given in Table III , together with those obtained in the model dependent analysis of 
Section IV . One can see that the values are similar and both the model independent and the model dependent analyses 
favor large mixing (0.1 < tan^i? < 10). This indication is very strong in both the model independent and the model 
dependent Global Analyses. 

Small mixing is obtained either through small (< 0.1) or large (> 10) values of t&v?d, that correspond, respectively, 
to i^e — vi and Ve — V2 (see Eq. (|l|)). Table III shows that both small and large values of tan^i? are disfavored, especially 
large values, in what could be called "extreme dark side" (tan^i? > 10), following the "dark side" denomination of the 
region tan^i? > 1 0. 
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Figures |l^ and show the posterior distribution of Am^ in the Rates Analysis and Global Analysis, respectively. 
The high and wide peak at Am^ ~ 10~^ — 10^^ eV^ in both figures implies that large values of Ato^ are favored. This 
is confirmed by the values of the integral probabilities of large (> 10~^eV^) and small (< 10^^ eV^) values of Am^ 



listed in Table [II, that are comparable to those obtained in the model dependent analysis, also reported in Table III. 

The model independent indications in favor of large mixing and large values of Am^ are very important for future 
terrestrial experiments that could explore this region of the parameter space with reactor neutrinos ifzsl, [79|. 



C. Neutrino Fluxes 



In this section we present the information on neutrino fluxes obtained through the marginal posterior distributions 



p($j|D,I) = J ]J d$,- J dC J dlog(tan2i?)dlog(ATO2)p(tan2i9, Am^$,C|D, 



I) 



(34) 



As could have been expected, we found that only the posterior distributions of the pp, ^Be and *B fluxes, to which 
the existing experiments are most sensitive, are significantly different from the corresponding prior distributions. 
The posterior distributions of all the other fluxes, including hep, is practically flat in the same range as the prior 
distribution. 

These distributions of the pp, ^Be and '^B fluxes are shown in Figs. [H- 



20 , In all these figures the thick solid line is 



the model independent posterior distribution, the thick dashed line is the prior distribution and the thick dotted line 
is the Gaussian BP2000 Standard Solar Model distribution. The intervals in which the thick solid line lies above the 
thin horizontal dotted, dashed and solid lines have, respectively, 90%, 95% and 99% probability. The thick solid lines 
representing the posterior distributions in Figs. |l5-20 have small fluctuations because the integrals over the fluxes 
and cross sections in Eq. (Q) have been performed through a Monte Carlo. The fluctuations of the thick dashed lines 
is due to the fact that the prior distribution is generated by the Monte Garlo. 

Table VII gives the credible intervals of the pp, ^Be and *B fluxes with 90%, 95% and 99% probability in the Rates 



and Global Analyses. 

Figures ^ and ^ show that the posterior distribution of the pp flux is in good agreement with the SSM distribution, 
with a slight prefe renc e for high values. The upper limit of the credible intervals with 90%, 95% and 99% probability 
reported in Table VII does not give a significant information, since it almost coincides with the upper limit of the 
prior distribution, 
hand, 




as explained in Section 



On the other 



that is determined by the luminosity constraint 
the lower limit of the credible intervals reported in Table [VII| is interesting, because it is significantly larger 
than the lower limit of the prior distribution, that follows from the prior distributions of the other fluxes and the 
luminosity constraint (|30|). In other words, the existing solar neutrino data allow to improve our knowledge of the pp 
flux, increasing its allowed lower limit. 

Also the posterior distributions of the ^Be flux shown in Figures |l^ and |l8| ar e in good agreement with the SSM 
distribution, with a slight preference for low values. One can see from Table VII that, according to the model 
independent analysis of solar neutrino data, the ^Be flux could be very small, but it is limited from above. Especially 
the Global Analysis restricts the allowed upper valu e for the ^Be flux significantly below the upper value of the prior 
distribution, as one can see from Fig. |l^ and Table VII. Let us remember that the ^Be neutrino flux has also been 
constrained in Re f. ISOl within ±25% (la error) of the SSM value from helioseismological measurements. 

Finally, Figures |l9| and 20 show that the posterior distribution of the ^B flux is in excellent agreement with the SSM 
distribution and it is bounded both from above and from below in an interval smaller than the prior distribution. This 
is a very interesting result, that agrees with the ranges for the *B flux obtained in Refs jS Jl g g |7[ |20[ |24[ H, H 
from the comparison of the SNO and Super-Kamiokande total rates only, and in Ref. p4|, p9| from the least-squares 
analysis of all solar neutrino data with the ^B flux considered as a free parameter. For example, the authors of Ref. 
found 0.44 < $SB/$fi'^ < 1.62 at 3ct, i.e. 99.73% confidence level, 
with 99.73% probability is 



We find that the credible interval of $83 /^gg 



0.42 < $SB/$f™ < 1.75 



in the Rates Analysis and 



0.45 < $83/$!™ < 1.42 



(35) 



(36) 



in the Global Analysis. Although the meaning of confidence level intervals and credible intervals is different, it is fair 
to say that the Global Analysis allows to restrict the allowed range of the *B neutrino flux more than the comparison 
of the SNO and Super-Kamiokande total rates only. 
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D. Oscillations versus no oscillations 

Bayesian Probability Theory allows to compare the compatibility of different hypotheses with the observed data 
calculating the ratio of their probabilities, that is usually called "odds" . Hence, we can compare in a model independent 
way the oscillation hypothesis with the no oscillation hypothesis, i.e. the Standard Model of neutrino physics. 

In order to perform this task we must reconsider Bayes Theorem in the form (^). The posterior probability of the 
oscillation hypothesis is given by 

p(OSC|D,I)= y"dlog(tan2i^)dlog(Am2) j J dC p{tan^^, Am^ ,<i>,C\J},l) 
_ /dlog(tan2i9)dlog(Am2) /d$ /dCp(D| tan^i?, Am^ $, C, I) p(tan2?9, Am^ji) p($|i) p(c|I) 



p(D|I) 

The flat prior distribution in the log(tan^i?)-log(Am^) plane of Figs. is given by 



(37) 



p(tan2^, Am^ |I) = ,^ , .. , (38) 

Alog(tan^w) Alog(ATO^) 

where Alog(tan^z9) and Alog(Am^) are the logarithmic ranges of tan^i? and Am^: Alog(tan^i?) = 7 and Alog(ATO^) = 
11. 

We compare the probability of neutrino oscillations (|37| ) with the probability of no oscillations 

fd$ rdCp(D|tan2|? = 0, Am2 = 0,$,C,I)p($|I)p(C|I) , , 

p(NO - OSC|D, I) = ^ ^ i-—^ ^-^^^1^^ ' ' ^ JFK \ li^y \ J _ ^3g^ 

Notice that in the ratio of Eqs. ( ^9| ) and ( |37| ) the unknown probability p(D|I) cancels. 
Our result for the comparison of no neutrino oscillations versus neutrino oscillations is 



p(NO-OSC|D,I) 
p(OSC|D,I) 

in the Rates Analysis and 

p(NO-QSC|D,I) 
p(OSC|D,I) 



8.7 X 10"'' (40) 



= 7.4x10"". (41) 



in the Global Analysis. Therefore, there is a rather strong model independent indication in favor of neutrino oscilla- 
tions, in agreement with those obtained in Refs. 0, ^ ^ ^, ^. The odds in favor of neutrino oscillations are slightly 
stronger in the Rates Analysis than in the Global^nalysis because the Super-Kamiokande energy spectrum is rather 
flat, which would be in agreement with the hypothesis of no oscillations if the Super-Kamiokande rate were equal to 
the SNO rate. 

Finally, let us notice that the prior distribution (|3^) and the result of the comparison of the oscillation and no- 
oscillation hypotheses that we have calculated depends on the ranges Alog(tan^i9) and Alog(ATO^). However, one 
must always keep in mind that these ranges have not been chosen randomly, but follow from our prior knowledge that 
there can be a measurable oscillation effect in solar neutrino experiments only if the neutrino mixing parameters lie 
within them. It is possible that somebody may argue in favor of a slight expansion or contraction of the prior ranges 
of Alog(tan^'i9) and Alog(ATO^) that we have considered. However, the effect of such expansion or contraction on the 
result of the comparison of the oscillation and no-oscillation hypotheses is small, not changing its order of magnitude. 

VI. CONCLUSIONS 

In this paper we have analyzed the data of solar neutrino experiment in terms of Ve J^/i,r oscillations without 
assuming the Standard Solar Model values of the neutrino fluxes. Working in the framework of Bayesian Probability 
Theory, the assumed lack of knowledge of the neutrino fluxes has been quantified through a flat prior probability 
distribution for the neutrino fluxes, that have been considered as unknown quantities to be determined, if possible, 
from the analysis of the data. 



In Section IV we have checked the reliability of our approach using a multi-normal prior probability distribution 



for the neutrino fluxes corresponding to the prediction of the BP2000 Standard Solar Model. We have shown that the 



14 



resulting allowed regions for the neutrino mixing parameters tan^i?, Am^ are practically equal with those obtained 
with the standard method in which the uncertainties of the fluxes are taken into account in the covariance matrix 
that determines the likelihood function. 

The results of the model independent analysis presented in Section ^ show that the present solar neutrino data allow 
to derive rather stringent model independent information on the neutrino mixing parameters, improving dramatically 
the conclusions obtained several years ago in another model independent analysis of solar neutrino data. 

Figures ^ and ^, obtained with the Rates Analysis and Global Analysis, respectively, show that the allowed regions 
of the neutrino mixing parameters are rather restricted and not much larger than those obtained in the SSM model 
dependent analysis (shown in Figs, ^and H). One can distinguish clearly the usual SMA, LMA, LOW and VO regions. 
Therefore, we conclude that the indication that the true values of the neutrino mixing parameters lie in one of these 
regions is robust. 

We have also shown in Section V B that there is a strong model independent indication in favor of large neutrino 
mixing and large values of Am^ . This is very important for future terrestrial experiments that could explore this part 



of the parameter space of neutrino mixing [ f78| , 79 . In particular, the model independent Global Analysis strongly 
disfavors the SMA solution of the solar neutrino problem. This is an impressive consequence of the high quality of 
the existing data. 

In Section V C we have presented the posterior distributions of the pp, ^Be and *B neutrino fluxes, to which the 
existing data are most sensitive (the posterior distribution of the other fluxes is similar to the prior distribution). 
We have shown that the posterior distributions of these fluxes are in excellent agreement with the BP2000 Standard 
Solar Model distribution, with a slight preference of a high pp flux and a low ''Be flux. The *B is severely constrained 
in an interval around the SSM value (see Table VII and Eqs. (^Sf), (|3^)), in agreement with the limits found in 
Refs. [||,|,|,g,|pl. 

Finally, in Section VD| we have shown that the hypothesis of no neutrino oscillations is strongly disfavored with 
respect to the hypothesis of i^e t'/i.r oscillations, in agreement with the results obtained in Refs. [||, ^, ^, ||, |^. 

In conclusion, we would like to emphasize the usefulness of the quantity and quality of the existing solar neutrino 
data, that nowadays allow to obtain stringent model independent information on neutrino physics and solar physics. 
We enthusiastically look forward towards the realization of new more sensitive solar neutrino experiments in order to 
improve our knowledge. 
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j 


Detection Material and Process 


Data 


1 


(Homestake |p|]) 


2.56 ±0.23 SNU 


2 


"Ga: !/e+ "Ga^"Ge + e" 
(GALLiiX p + GNO ||1(| + sAGli ||!i|J) 


74.7 ± 5.1 SNU 


3 


D2O: Ue+d-^p + p + e^ 
(SNO 1) 


0.347 ± 0.028 


4 


H2O: + e ^ + e 
(Super-Kamiokande [|l)) 


0.459 ± 0.017 



TABLE I: The rates measured in solar neutrino experiments. The rates of the Homestake and GALLEX+SAGE+GNO 
experiments are expressed in SNU units (1 SNU = 10~^^ events atom~^ whereas the results of the Kamiokande and SNO 
experiments are expressed in terms of the ratio of the experimental rate and the BP2000 Standard Solar Model prediction [ ^rf . 
The statistical and systematic uncertainties have been added in quadrature. The GALLEX+SAGE+GNO rate is a weighted 
average of the GALLEX+GNO rate reported in Ref. @ and the SAGE rate reported in Ref. ^. The rate of the SNO 
experiment is that measured through CC weak interactions. 



b 
day 


b 

day 


^3b 

(cm^ s) 


(cm^ s) 


D.cor 


fa, cor 


1 


20 


6.21 X 10"" 


1.98 X 10" ' 


0.0039 


0.0020 


2 


21 


6.62 X lO"** 


1.98 X 10" ' 


0.0039 


0.0020 


3 


22 


7.09 X 10"** 


1.98 X 10" ' 


0.0047 


0.0029 


4 


23 


7.63 X 10"'' 


1.98 X 10" ' 


0.0055 


0.0058 


5 


24 


8.22 X lO"** 


1.98 X 10" ' 


0.0062 


0.0077 


6 


25 


9.23 X 10"" 


1.98 X 10" ' 


0.0070 


0.0106 


7 


26 


1.00 X 10"' 


1.98 X 10" ' 


0.0078 


0.0134 


8 


27 


1.12 X lO"' 


1.98 X 10" ' 


0.0086 


0.0163 


9 


28 


1.29 X 10"' 


1.98 X 10" ' 


0.0093 


0.0203 


10 


29 


1.47 X 10"' 


1.98 X 10" ' 


0.0101 


0.0242 


11 


30 


1.72 X 10"' 


1.98 X 10" ' 


0.0109 


0.0292 


12 


31 


2.03 X 10"' 


1.98 X 10" ' 


0.0117 


0.0331 


13 


32 


2.44 X 10"' 


1.98 X 10" ' 


0.0124 


0.0381 


14 


33 


2.91 X 10"' 


1.97 X 10" ' 


0.0132 


0.0440 


15 


34 


3.72 X 10"' 


1.97 X 10" ' 


0.0140 


0.0500 


16 


35 


4.72 X 10"' 


1.97 X 10" ' 


0.0148 


0.0570 


17 


36 


6.08 X 10"' 


1.97 X 10" ' 


0.0155 


0.0640 


18 


37 


8.21 X 10"' 


1.97 X 10" ' 


0.0163 


0.0730 


19 


38 


1.95 X 10"" 


1.94 X 10" ' 


0.0205 


0.1058 



TABLE II: Values of the relative theoretical uncertainties 5, 
energy bins used in the Global Analysis (see Section II B). 
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and experimental uncertainties S^^^^l of the Super-Kamiokande 
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BP2000 SSM 


Model Independent 


Rates 


Global 


Rates 


Global 


SMA 


9.9 X 10"^ 


7.2 X 10"^ 


0.15 


1.0 X 10"^ 


LMA 


0.72 


0.86 


0.67 


0.80 


LOW 


6.1 X 10"^ 


0.12 


3.8 X 10"^ 


0.13 


VO 


0.12 


1.9 X 10"^ 


0.14 


6.6 X 10"^ 


tan^i? < 0.1 


9.9 X 10"^ 


7.2 X lO"'^ 


0.18 


8.2 X 10"=* 


0.1 < tan^'d < 10 


0.90 


1.00 


0.82 


0.99 


tan^j? > 10 


6.5 X 10"* 


1.3 X 10"* 


1.1 X 10"^ 


1.7 X 10"^ 


Am^ < 10~^eV^ 


0.32 


0.14 


0.41 


0.20 


Am^ > 10~^eV^ 


0.68 


0.86 


0.59 


0.80 



TABLE IIL Integral probabilities of regions in the tan^i?-Am^ plane and of intervals of tan^i? and Am'^ obtained with the 
model dependent and model independent Rates Analysis and Global Analysis. The SMA, LMA, LOW and VO regions are 
defined in Eqs. 



i 


L 

PP 


2 

pep 


3 
hep 


4 

•^Be 


5 

8B 


6 

13n 


7 

15o 


8 




+0.094 
-0.096 


±1 


+9 
-1 


±1 


±1 


±1 


±1 


±1 




iO.OLO 


±0.015 




±0.100 


±0.200 


±0.210 


±0.250 


±0.250 


NACRE (3.842) [ 


Ll 




0.002 


-0.007 


0.011 


0.008 


0.077 


-0.111 


-0.129 


-0.059 


BP2000 (3.844) [ 


LI 




0.002 


0.000 


0.000 


0.010 


0.020 


0.015 


0.017 


0.018 


NACRE (3.844) [ 


Ll 




0.003 


-0.007 


0.011 


0.017 


0.097 


-0.100 


-0.117 


-0.043 


ASOO [ 


11 




0.007 


0.007 


0.011 


-0.031 


-0.069 


-0.042 


-0.050 


-0.053 


GN 93 


11 




-0.002 


-0.007 


-0.011 


0.023 


0.051 


0.128 


0.135 


0.155 


Pre-M.S 


|1 




0.000 


-0.007 


-0.011 


0.021 


0.048 


0.124 


0.131 


0.149 


Rotation [] 


i 


0.005 


0.000 


-0.011 


-0.019 


-0.02S 


0.016 


0.015 


0.028 


Radius78 [: 


i 


-0.002 


-0.007 


-0.011 


0.023 


0.051 


0.128 


0.135 


0.155 


Radiussos 


Lil 


-0.002 


-0.007 


-0.011 


0.023 


0.051 


0.128 


0.135 


0.155 


No Diffusion 




1 


0.0L7 


0.021 


0.032 


-0.117 


-0.234 


-0.254 


-0.279 


-0.2S1 


Old Physics 


\\1 


1 


0.000 


0.007 


-0.011 


0.029 


0.020 


0.053 


0.048 


0.052 


S3,4 = |l 


i) 


0.076 


0.107 


0.086 


-1.000 


-1.000 


O.lSl 


0.175 


0.190 


Mixed 






0.030 


-0.093 


-0.333 


-0.252 


-0.182 


-0.445 


-0.365 


-0.359 


BP98 1 






-0.002 


-0.007 


-0.774 


0.006 


0.020 


0.104 


0.108 


0.124 


BP95 1 


82 




-0.007 


0.000 


-0.871 


O.OSO 


0.311 


0.128 


0.135 


0.151 


DS94 


0.017 


0.039 


-0.866 


-0.107 


-0.29S 


-0.074 


-0.0S3 


-0.211 


TL93 |3 


, \ 


J 


0.013 


-0.007 




-0.090 


-0.123 


-0.301 


-0.337 




BP92 1 


85 




0.008 


0.021 


-0.86S 


0.025 


0.127 


-0.102 


-0.113 


-0.043 


BUSS 1 


86 




0.008 


0.000 


-0.185 


-0.017 


0.141 


0.111 


0.0S8 


-0.0S3 



TABLE IV: Relative uncertainties M>i/^r^ of the prior distributions of the neutrino fluxes in the model independent analysis 
confronted with the relative uncertainties (A"l>f ^^/$f in the BP2000 Standard Solar Model and the differences between 
the fluxes predicted by several models and the BP2000 SSM. 
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Rates Analysis 



tan^i? 


Am^ (eV^) 


PP 


pep 


hep 


^Be 


«B 


13 N 


15q 




1.7 X 10"^ 


8.0 X lO^'* 


0.935 


0.220 


8.907 


1.992 


1.290 


1.090 


0.177 


0.063 


0.20 


1.2 X 10~'^ 


0.962 


0.960 


6.752 


1.508 


1.354 


1.362 


0.764 


1.253 


0.49 


1.4 X lO"'^ 


0.970 


0.813 


1.012 


1.437 


1.003 


0.805 


0.749 


1.250 


3.59 


8.0 X 10"" 


1.031 


1.779 


1.449 


0.739 


0.904 


0.800 


0.033 


1.032 


1.0 X 10"^ 


3.0 X 10"^ 


1.035 


1.344 


5.116 


0.349 


1.981 


1.967 


1.861 


1.632 


5.0 X 10^2 


1.0 X 10"'' 


0.916 


1.885 


5.203 


1.996 


1.997 


1.758 


1.173 


1.844 


1.00 


9.0 X 10"" 


0.965 


0.841 


3.686 


1.464 


1.563 


0.454 


1.015 


0.202 


1.00 


5.5 X 10"^^ 


0.983 


0.757 


0.040 


1.212 


0.446 


1.822 


0.853 


0.025 








1.089 


0.679 


0.249 


0.021 


0.428 


0.233 


0.103 


0.481 



TABLE V: Values of the neutrino fluxes normalized to the BP2000 SSM fluxes corresponding to the lines in Fig. ^ obtained 
in the Rates Analysis with the nine selected values of tan^i? and Am'^ listed in the first two columns. The first four rows 
correspond, respectively, to the best-fit points in the SMA, LMA, LOW and VO regions in the Rates Analysis. The fifth to 
eighth rows correspond to points in excluded areas in Fig. p^. The last row corresponds to the case of no oscillations. 



Global Analysis 



tan^i? 


Am^ (eV") 


PP 


pep 


hep 


^Be 


«B 


13 N 


15Q 


ITp 


1.7 X 10"^ 


8.0 X 10"*' 


0.929 


0.165 


3.922 


1.862 


1.260 


0.262 


1.523 


1.557 


0.20 


1.2 X 10"^ 


0.966 


1.414 


0.335 


1.507 


1.286 


1.906 


0.301 


0.699 


0.49 


1.4 X lO"'^ 


0.968 


0.045 


1.561 


1.520 


0.992 


0.617 


0.503 


1.147 


3.59 


8.0 X 10~" 


1.014 


0.987 


1.004 


0.896 


0.931 


0.048 


0.651 


0.426 


3.8 X 10"'' 


5.3 X 10"^ 


1.035 


0.697 


2.764 


0.535 


0.526 


0.394 


1.186 


0.114 


0.34 


6.1 X lO"-'' 


1.004 


1.248 


5.544 


1.047 


1.104 


1.854 


0.186 


0.433 


0.60 


1.5 X 10"'^ 


0.994 


0.147 


5.910 


1.223 


0.915 


0.319 


0.342 


0.967 


3.16 


5.0 X 10"'° 


1.050 


1.111 


4.100 


0.510 


0.620 


0.008 


0.141 


0.255 








1.088 


0.065 


0.664 


0.071 


0.438 


0.032 


0.031 


0.340 



TABLE VI: Values of the neutrino fiuxes normalized to the BP2000 SSM fiuxes corresponding to the lines in Fig. obtained 
in the Global Analysis with the nine selected values of tan^i? and Am'^ listed in the first two columns. The first four rows 
correspond, respectively, to the best-fit points in the SMA, LMA, LOW and VO regions in the Rates Analysis. The fifth to 
eighth rows correspond, respectively, to the best-fit points in the SMA, LMA, LOW and VO regions in the Global Analysis. 
The last row corresponds to the case of no oscillations. 
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Rates Analysis 


Global Analysis 


90% 


95% 


99% 


90% 


95% 


99% 




[0.95, 1.08] 


[0.93, 1.08] 


[0.92, 1.09] 


[0.99, 1.08] 


[0.98, 1.08] 


[0.95, 1.09] 




[0.02, 1.65] 


[0.02, 1.82] 


[0.02, 1.95] 


[0.02, 1.15] 


[0.02, 1.28] 


[0.02, 1.62] 




[0.62, 1.48] 


[0.58, 1.58] 


[0.48, 1.68] 


[0.62, 1.22] 


[0.58, 1.28] 


[0.48, 1.35] 



TABLE VII: Credible intervals for the pp, "^Be and fluxes with 90%, 95% and 99% probability obtained with the model 
independent Rates and Global Analyses. 
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tan^ d 



FIG. 1: Credible regions obtained with the Rates Analysis of Homestake, GALLEX+GNO+SAGE, SNO and Super- 
Kamiokande total rates, assuming the BP2000 Standard Solar Model neutrino fluxes. 
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FIG. 2: Credible regions obtained with the Global Analysis of Homcstako, GALLEX+GNO+SAGE and SNO rates and 
Super-Kamiokande day and night spectra, assuming the BP2000 Standard Solar Model neutrino fluxes. 
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FIG. 3: Credible regions obtained with the model independent Rates Analysis of Homestake, GALLEX+GNO+SAGE, SNO 
and Super-Kamiokande total rates. 
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FIG. 4: Credible regions obtained with the model independent Rates Analysis of Homestake, SNO and Super-Kamiokande 
data, neglecting the Gallium data of the GALLEX, GNO and SAGE experiment. 
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tan^ d 



FIG. 5: Credible regions obtained with the model independent Rates Analysis of GALLEX+GNO+SAGE, SNO and Super- 
Kamiokande data, neglecting the Chlorine data of the Homestake experiment. 




FIG. 6: Credible regions obtained with the model independent Rates Analysis of Homestake, GALLEX+GNO+SAGE and 
Super-Kamiokande data, neglecting the rate measured in the SNO experiment. 




FIG. 7: Credible regions obtained with the model independent Rates Analysis of Homestake, GALLEX+GNO+SAGE and 
SNO data, neglecting the rate measured in the Super-Kamiokande experiment. 
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FIG. 8: Experimental rates (points with errorbars) and theoretical predictions for nine selected values of tan^i? and Am'^ in 
the Rates Analysis. (A). Solid line: tan^i? = 1.7 x 10"^ and Am^ = 8.0 x 10"^ eV^ best fit point in the SMA region in the 
Rates Analysis; long-dashed line: ta.v^'d = 0.20 and Am^ = 1.2 x 10~^eV^, best fit point in the LMA region in the Rates 
Analysis; short-dashed line: tax^d = 0.49 and Am^ = 1.4 x 10~'^eV^, best fit point in the LOW region in the Rates Analysis; 
dotted line: tan^i? = 3.59 and Am^ = 8.0 x 10~^^ eV^, best fit point in the VO region in the Rates Analysis. (B). Solid line: 
tan^i? = 1.0 X 10"^ and Am^ = 3.0 x 10"^ eV^; long-dashed line: tan^-d = 5.0 x 10"^ and Am^ = 1.0 x 10"** cV^; short-dashed 
line: tan^i? = 1.00 and Am^ = 9.0 x 10"" eV^; dotted line: tan^i? = 1.00 and Am^ = 5.5 x 10"^^ eV^; dash-dotted line: no 
oscillations. 
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FIG. 9: Credible regions obtained with the model independent Global Analysis of Homestake, GALLEX+GNO+SAGE and 
SNO rates and Super-Kamiokande day and night spectra. 
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FIG. 10: Experimental rates (points with crrorbars) and theoretical predictions for nine selected values of tan^i? and Am^ 
in the Global Analysis. The index j labels the solar data points: j = 1,2,3 for the Gallium, Chlorine and SNO event rates, 
respectively; j = 4, . . . , 22 for the Super-Kamiokande energy spectrum in the day; j = 23, . . . , 41 for the Super-Kamiokamde 
energy spectrum in the night. (A). Solid Une: tan^?? = 1.7 x 10~^ and Am^ = 8.0 x 10"'^ eV^, best fit point in the SMA region 
in the Rates Analysis; long-dashed line: tan^t? = 0.20 and Am^ = 1.2 x 10~*eV^, best fit point in the LMA region in the 
Rates Analysis; short-dashed line: ta.v?{) = 0.49 and Am'^ = 1.4 x 10^^ cV^, best fit point in the LOW region in the Rates 
Analysis; dotted lino: tan^i? — 3.59 and Am'^ — 8.0 x 10^^^ eV'^, best fit point in the VO region in the Rates Analysis. (B). 
Solid line: tan^j? = 3.8 x 10""* and Am^ = 5.3 x 10"® eV^ best fit point in the SMA region in the Global Analysis; long-dashed 
line: tan^i? = 0.34 and Am^ = 6.1 x 10~^eV^, best fit point in the LMA region in the Global Analysis; short-dashed line: 
tan^?? = 0.60 and Am^ = 1.5 x 10~'^eV^, best fit point in the LOW region in the Global Analysis; dotted Une: tan^?? = 3.16 
and Am^ = 5.0 x 10~^° eV^, best fit point in the VO region in the Global Analysis; dash-dotted line: no oscillations. 




FIG. 11: Posterior distribution of tan^i) obtained with the model independent Rates Analysis. The intervals in which the thick 
solid curve lies above the horizontal dotted, dashed, and solid lines have, respectively, 90%, 95% and 99% posterior probability 
to contain the true value of tan^i?. 




FIG. 12: Posterior distribution of tan^i? obtained with the model independent Global Analysis. The intervals in which the 
thick solid curve lies above the horizontal dotted, dashed, and solid lines have, respectively, 90%, 95% and 99% posterior 
probability to contain the true value of tan^i?. 
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FIG. 13: Posterior distribution of Arn^ obtained with the model independent Rates Analysis. The intervals in which the thick 
solid curve lies above the horizontal dotted, dashed, and solid lines have, respectively, 90%, 95% and 99% posterior probability 
to contain the true value of Am^. 




FIG. 14: Posterior distribution of Am'^ obtained with the model independent Global Analysis. The intervals in which the thick 
solid curve lies above the horizontal dotted, dashed, and solid lines have, respectively, 90%, 95% and 99% posterior probability 
to contain the true value of Am^. 
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FIG. 15: The thick solid line is the posterior distribution of the pp flux obtained with the model independent Rates Analysis. 
The intervals in which the thick solid line lies above the horizontal thin dotted, dashed, and solid lines have, respectively, 90%, 
95% and 99% posterior probability to contain the true value of the pp flux. The gaussian thick dotted curve is the BP2000 
SSM distribution and the thick dashed curve is the prior distribution in the model independent analysis. 



36 




FIG. 16: The thick solid hne is the posterior distribution of the pp flux obtained with the model independent Global Analysis. 
The intervals in which the thick solid line lies above the horizontal thin dotted, dashed, and solid lines have, respectively, 90%, 
95% and 99% posterior probability to contain the true value of the pp flux. The gaussian thick dotted curve is the BP2000 
SSM distribution and the thick dashed curve is the prior distribution in the model independent analysis. 
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FIG. 17: The thick solid line is the posterior distribution of the '^Be flux obtained with the model independent Rates Analysis. 
The intervals in which the thick solid line lies above the horizontal thin dotted, dashed, and solid lines have, respectively, 90%, 
95% and 99% posterior probability to contain the true value of the pp flux. The gaussian thick dotted curve is the BP2000 
SSM distribution and the thick dashed curve is the prior distribution in the model independent analysis. 
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FIG. 18: The thick solid line is the posterior distribution of the '^Be flux obtained with the model independent Global Analysis. 
The intervals in which the thick solid line lies above the horizontal thin dotted, dashed, and solid lines have, respectively, 90%, 
95% and 99% posterior probability to contain the true value of the pp flux. The gaussian thick dotted curve is the BP2000 
SSM distribution and the thick dashed curve is the prior distribution in the model independent analysis. 
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FIG. 19: The thick sohd line is the posterior distribution of the flux obtained with the model independent Rates Analysis. 
The intervals in which the thick solid line lies above the horizontal thin dotted, dashed, and solid lines have, respectively, 90%, 
95% and 99% posterior probability to contain the true value of the pp flux. The gaussian thick dotted curve is the BP2000 
SSM distribution and the thick dashed curve is the prior distribution in the model independent analysis. 
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FIG. 20: The thick solid line is the posterior distribution of the flux obtained with the model independent Global Analysis. 
The intervals in which the thick solid line lies above the horizontal thin dotted, dashed, and solid lines have, respectively, 90%, 
95% and 99% posterior probability to contain the true value of the pp flux. The gaussian thick dotted curve is the BP2000 
SSM distribution and the thick dashed curve is the prior distribution in the model independent analysis. 



